Non-markovian global persistence in phase-ordering 
kinetics 

Malte Henkel 1 and Michel Pleimling 2 

1 Groupe de Physique Statistique, 

Departement de Physique de la Matiere et des Materiaux, Institut Jean LamouiQ, 
CNRS - Nancy Universite - UPVM, B.P. 70239, 
F - 54506 Vandoeuvre les Nancy Cedex, France 

2 Department of Physics, Virginia Polytechnic Institute and State University, 
Blacksburg, Virginia 24061-0435, USA 

E-mail: henkel@lpm.u-nancy . fr, Michel . PleimlingOvt . edu 

Abstract. The persistence probability P s (t) of the global order-parameter of a 
simple ferromagnet undergoing phase-ordering kinetics after a quench from a fully 
disordered state to below the critical temperature, T < T c , is analysed. It is argued 
that the persistence probability decays algebraically with time in the entire low- 
temperature phase. For Markov processes, the associated global persistence exponent 
g = [2Xc — a\/(2z) is related to the autocorrelation exponent Ac- This relationship 
is confirmed for phase-ordering in the exactly solved II? Ising model and the d- 
dimensional spherical model. For the 2D Glauber-Ising model, the temperature- 
independent estimate 9 g = 0.063(2) indicates that the dynamics of the global order- 
parameter is described by a non-Markovian process. 
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Persistence properties of a stochastic process have been studied in a variety of 
dynamical systems, ranging from magnetic systems [U [2J, [31 SI El El [Tl [HI [9] to soap 
bubbles [10], from reaction-diffusion systems pj] to fluctuating interfaces [12], turbulent 
liquid crystals [13] and non-equilibrium surface growth [H] . The persistence probability 
Pit) is defined in a very general way as the probability that a given stochastic variable 
(observable) retains a characteristic feature over the time t. One commonly studied 
problem involves the question whether the stochastic variable, which can be a local or 
a global quantity, maintains its sign up to the time t. 

Over the years, the persistence properties of the fluctuations around the mean value 
of the local order parameter <p(t) have attracted much interest. Problems of this kind 
were initially formulated for the probability Pe{t) that the local order-parameter <p(t, r) 
has not yet changed its sign up to time t. For large enough times, one expects a power- 
law decay Pe(t) ~ t~ 9e if the system is at an equilibrium critical point, and where 0£ is 
the local persistence exponent. This has been studied for simple diffusion [151 [16] an d for 
one-dimensional exactly solvable models quenched to T = pD, El HI El [6] , for example 
in the g-states Potts model one has the remarkable result 0t — — | + [| arccos(-^=p] 2 /2 
[5]. Values of 0g for higher-dimensional models have been extracted from numerical 
simulations, see [3] for data on the 2D Ising model. 

It has turned out to be very useful [7J to study the persistence properties of the 
global order-parameter 



averaged over spatial domains Q with large volumes 1 <C <C \V\, where |V| is 
the total volume of the system. The global persistence probability P s (t) is defined as 
the probability that 4>o(t) has not changed sign up to time t. For large times, one 
expects a power-law decay P g (t) ~ t~ dg which defines the global persistence exponent 
9g. This scenario was shown to apply in particular to non-equilibrium critical dynamics 
(T = T c ) [7J. Here, an application of the central limit theorem shows that 0o(^) 
should be described by a Gaussian stochastic process and if that process is furthermore 
Markovian, one has the scaling relation 9 g z = Xq — d+1 — r]/2 [7J. The definition 
of the exponents 77 and Xc will be recalled below. For critical systems, this relation is 
satisfied in a few exactly solvable systems such as the ID Glauber-Ising model quenched 
to T = or the spherical model in 2 < d < 4 dimensions (for all of which the dynamical 
exponent happens to be z = 2). However, the scaling relation does break down in more 
generic systems, where z ^ 2 in all known examples. One concludes [7J that the process 
describing the long-time collective dynamics of <po{t) and which arises from renormalised 
field-theory, is in general non-Markovian, see [IT] for a review. This has been confirmed 
by a large variety of analytical and simulational studies, either for critical ferromagnets 
(including effects of disorder) [18] or else for genuine non-equilibrium systems at an 
absorbing phase transition [19]. More recent developments consider the cross-over in 
the global persistence as a function of the initial value of the magnetisation [20] or the 
behaviour in the vicinity of surfaces [2H [22] or on sub-manifolds [23] . 




(1) 
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Aging phenomena, which are observed in a broad variety of systems with slow 
relaxation dynamics [2H [25], (26J [27] 12H] , are paradigmatic examples of systems that are 
far from stationarity, and of great potential for applications. Initially found in glassy 
systems, they also occur in simple ferromagnets, which are often considerably easier 
to analyse. This kind of non-equilibrium collective phenomena may be conveniently 
realized by quenching a system from a fully disordered initial state to a temperature 
T < T c , where T c is the critical temperature. The specific aspects of aging, namely (i) 
slow, non-exponential relaxation, (ii) breaking of time-translation-invariance and (iii) 
dynamical scaling, are conveniently studied through the behaviour of two-time quantities 
such as the two-time correlator C(t, s;r) = (<ft(t, r)0(s, 0)) of the order parameter. 
For simple ferromagnets, one finds a dynamical scaling regime when both times are 
sufficiently large and if t — s ^> t m icro, where t m i cr o is a microscopic reference time. Then 

C(t,s;v) = S - b f c (t/s,\v\ z /(t-s)) (2) 

where 6 = for T < T c and b = (d — 2 + r])/z for T = T c , where r] describes the decay of 
the critical equilibrium correlations and z is the dynamical exponent. For y = t/s large, 
one has fc(y,0) ~ y~ x °/ z where Ac is the autocorrelation exponent [29] [30], whose 
value depends on whether T < T c or T = T c . 

In this work, we wish to return to an analysis of the properties of the global 
persistence for simple ferromagnets which undergo phase-ordering kinetics after a quench 
from a fully disordered initial state, with a vanishing average initial global magnetisation 
(0o(O)) = 0, to T < T c with T c > oJ§| It has been known since a long time that z = 2 
|31[ [32] and current theories for the growth laws in phase-ordering kinetics usually start 
from a local differential equation dt<fi = —S J 7 /8(f) where T is a local Landau- Ginzburg 
functional. Since the temperature T < T c is thought to be irrelevant, one sets T = and 
analyses the long-time behaviour of the solutions of the above equation, after having 
averaged over the totally disordered initial state [31], [32] [21] . Although this set-up is 
perfectly local in time, the behaviour of the global persistence probability is non-trivial 
in general. We find: 

(i) The global order-parameter 4>o(t) is described by a Gaussian stochastic process and 
its global persistence probability P g (t) ~ t~° 9 decays algebraically. 

(ii) If the process for <po{t) is also Markovian, then 9 g = #™ ark , where [8] [9] 

9™ A z = \ c -d/2>0 (3) 

The bound follows from the Yeung-Rao-Desai inequality [33J. 

(iii) Eq. ([3]) is satisfied for several exactly solvable systems such as the ID Glauber-Ising 
model at T = [7] and the spherical model at T < T c and for d > 2. These systems 
are explicitly Markovian. 

§ For an ordered initial state, however, one expects P g (t) — > Pg,oc(T) exponentially fast in t and 
without dynamical scaling. For T < T c , there is no analogue to the cross-over observed in [5D] for 
critical quenches. 
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(iv) While the results stated so far come from essentially standard procedures, we also 
performed simulations of the 2D Glauber-Ising model quenched to T < T c « 2.27 
and find 

f 0.062(2) ; for T = 1.0 
9 \ 0.065(2) ; for T = 1.5 ' 1 1 

These results do not agree with eq. ([3]). Hence, the process <po(t) is not Markovian 
in general. This example of a presumably non-integrable system also illustrates that 
the simple value z = 2 of the dynamic exponent, as it holds true in phase-ordering 
kinetics [311 32J, does not imply the Markov property. 

Generalising the lines of thought developed in [7J to quenches below T c [8j [5] , consider 
the formation of ordered domains. Shortly after the quench, they should be rather small, 
with a linear size £o at most of the order of the lattice constant. For later times, the 
firmly established dynamical scaling [21] (for a proof of dynamical scaling in 2D phase- 
ordering, see [31]) leads to a single relevant time-dependent length-scale L(t) ~ t l l z 
and the dynamics of the system is described by the movement of the domain walls of 
independent domains of linear sizes larger than L(t), where the volume |V| ^> L(t) d 
of the system must be large enough. For any given time t, the first two moments of 
the global order-parameter are finite and read (<j)o(t)) — and (0n(i)) ~ t d l z . Since 
one averages over the configuration space of disordered initial states the central limit 
theorem is applicable, hence the stochastic process describing the time-evolution of 4>o(t) 
should be Gaussian, also for T < T c [HIE]. Again because of dynamical scaling, the two- 
time global correlator is (4> (t)(j) (s)) = s d / z fc(t/s) with t > s and the asymptotic 
behaviour fc{y) ~ y^ x c)l z as y — > oo. Therefore the normalised global autocorrelator 
(t>s) 

N{t,s):= 0^)1 =f N (t/s) (5) 



has the asymptotic form fN(y) ~ y~( x c- d / 2 )/ z as y 3> 1- Changing temporal variables 
according to T = hit, the process describing </>o becomes stationary. According to 
Doob's lemma [36J, in this setting the Markov property is equivalent to having the 
exact correlator N(t, s) = N(T, S) = exp(— n\T — S\), with n = (2Ac — d)/{2z). Then 
Slepian's formula [37] shows that the global persistence probability indeed decays with 
time according to P g (t) ~ £- 6 s nark suc h that eq. ([3]) holds true. 

Not surprisingly, eq. ([3]) can be confirmed in some exactly solvable models the 
dynamics of which are explicitly Markovian. As a first example, we consider the ID 
Glauber-Ising model quenched to T = from a fully disordered state. Since T c = 
in ID, the global persistence in this model was already considered as an example for 
a critical quench [7J, but we shall see that the same result may also be interpreted in 
a low-temperature setting. Indeed, from the well-known results Ac* = 1 and z = 2 

|| If that average is not taken, simple scaling may be replaced by multi-scaling, see e.g. [3S] 
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eq. ([3]) gives the Markovian prediction #™ ark = A. The exact calculation gives the 
expected power-law decay and 9 g = \ [7], as it should be. 

As a second example, we consider the spherical model in d > 2 dimensions 
quenched to T < T c . The spherical model may be formulated in terms of real- valued 
spin variables S r subject to the constraint J2 r = with the total number J\f of 
lattice sites. The Hamiltonian describes ferromagnetic nearest-neighbour interactions 
and reads H = —J2(r,r')^rS T '- The non-conserved dynamics is given by a Langevin 
equation, with a totally disordered initial state. Following the usual lines of its 
exact solution [39], the two-time spin-spin correlator at momentum q reads Cq(t, s) = 



CJO) + 2T f dre Mq)T g{r) 
Jo 



Cq(s)e- UJ ^ t - s ^g(s)/g(t) where 

A M _ exp(-2cu(q)t) 

g{t) =f(t) + 2T fdrf{t-r)g{r) (6) 
Jo 

with fit) = (e~ 4 *J (4t)) d and Jo is a modified Bessel function, while uj(q) = 25Zf=i(l — 
cosgj) ~ q 2 . Then the normalised global correlator becomes 



N(t,s)- do ^ S 



yJC {t)6 (s) \ C (t)g(t) 



C (s)g(s) ^ 



For a fully disordered initial state C o (0) = 1. 

Indeed, for zero temperature T = 0, we readily have N(t,s) = 1, in agreement 
with the explicit Markov property in this model. In order to show that this result 
remains true for all < T < T c , we write first Co(t)g(t) = 1 + 2T /pdr g{r) and then 
solve the Volterra integral equation d^J) [391 HO] , which leads for large enough times to 
f dr g(r) = go + git l ~ d l 2 , where the explicitly known values of the constants #0,1 will 
not be needed. Hence 

N(t,s) = l + 0(s^ d/2 ) (8) 

indeed follows an exact power law in t/s, up to corrections to dynamical scaling, and 
Slepian's formula [37] gives 6 g = 0. Since the autocorrelation exponent = d/2 |39j, 
this agrees with the Markovian prediction (j3J), as expectedj^l 

Finally, we study the global persistence in the 2D Glauber-Ising model, quenched 
to T < T c from a fully random initial state with zero magnetisation. In Fig. [T] we 
show the global persistence probability P g (t), and also compare results for different 
values of T. The global persistence probability is simply the probability that the global 
magnetisation has not crossed zero until time t. Finite-size effects were monitored 
carefully by simulating systems of various system sizes. The data in Fig. [T] have been 
obtained for a system with 400x400 sites and result from an average over 80,000 different 

If A similar conclusion holds true for the long-range spherical model, with Hamiltonian TL = 
~Er#r' J ( r - *')SrS r > and J(r) - |r|- d - CT with < a < 2, quenched to T < T c . For non-conserved 
dynamics, one has z = a, Xc — d/2 [41] . and a straightforward extension of the above calculation shows 
that N(t, s) — 1, up to finite-time corrections, and that eq. ([3|) is satisfied. 
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Figure 1. Global persistence probability for the 2D Glauber-Ising model, quenched 
either to T — 1.0 or to T = 1.5. After an initial regime at short times a power-law 
regime is observed, with an exponent 9 g w 0.063 (dashed lines). Data were averaged 
over 80, 000 different runs with different initial states and different realizations of the 
noise. 

runs where we averaged over different initial conditions [35] and different realizations 
of the noise. Clearly, for sufficiently large times the decay of the global persistence is 
described by a power-law. The values of 9 g extracted are listed in eq. (pO). However, 
since = 1.25(2) [261 I2H], the estimates (j3J) are quite distinct from the Markovian 
prediction #™ ark = 0.125(2). Therefore, the effective long-time dynamics of the global 
order-parameter should not be described by a Markov process. This is our main result. 

The data in Fig. [T] seem to indicate the existence of an initial short-time regime, 
followed by the power-law decay. In fact, as we start from a disordered initial state, the 
dynamical correlation length L(t) increases as a function of time. One then expects that 
the power-law behaviour characteristic of the T = fixed point should dominate the 
persistence probability only once Lit) ^> £ where £ is the equilibrium correlation length 
which sufficiently close to T c should scale as £ ~ (T c — T)~ v . As long as L <C L{t) <C £, 
where L is a microscopic length scale, the leading behaviour of P g (t) should be the same 
as at the critical point. In order to verify this, we monitored the initial time behaviour as 
a function of temperature and found that the initial regime can be described by a power- 
law decay with a temperature-dependent effective exponent. The value of this exponent 
(as well as the cross-over time between the short-time regime and the late-time regime, 
see Fig. [1]) increases with increasing T. We find the values 0.18 for T = 1.8 and 0.20 
for T = 2.0, already close to the value 6 g = 0.237(3) (TJ |42j [20] of the global persistence 
exponent at the critical point. If L(t) is of the order of L , the evolution will depend 
on all sorts of non-universal, microscopic details, and no general, model-independent 
statement can be made. 

In Fig. [2] we show the normalised correlator N(t, s) for different waiting times s 
at the temperature T = 1.5. We first verify in Fig. [2^l that N(t,s) is indeed only a 
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Figure 2. (Colour online) Normalised correlator N(t, s) for the 2D Glauber-Ising 
model, quenched to T = 1.5. (a) Data obtained for different waiting times fall on 
a common master curve when plotted as a function of t/s. The effective exponent 
at large times is w 0.115, already close to the expected value 0.125. (b) Plotting 
(t/s) ' N(t,s) over against t/s shows that even for the largest accessible values of 
t/s, the asymptotic power-law regime is not yet completely reached. The error bars, 
which result from averaging over at least 5,000 independent runs, are much smaller 
than the symbol sizes. 

function of t/s, see fl5]). As we are exploring the time-dependence of the correlator up to 
t/s = 60, we have to simulate rather large systems with 800 x 800 spins in order to avoid 
finite-size effects. While in the accessed time-window the value of the effective exponent 
of the asymptotic decay ss 0.115 is already close to the expected value 0.125, we also find 
that this correlator does not follow a simple power law which gives additional evidence 
against O being described by a Markov process. 

Our result ([4]) for 6 g is consistent with data on the block persistence at T = 
[H E], where 9 g ~ 0.09 in the 2D Ising model and 9 g ~ 0.06 in the 2D time-dependent 
Ginzburg-Landau equation was found, in agreement with the expectation that T should 
be an irrelevant parameter. 

We observe that while for quenches to T = T c , all known systems give 6 g > #™ ark , 
see e.g. [3 UHl [TTl [201 EH E21 EH S2] , the presently available evidence suggests that the 
opposite inequality 9 g < #™ ark might hold true for quenches into the low-temperature 
phase T <T C (here, we exclude cases with T c = from the discussion). 

Summarising, we have studied the extension of a scaling relation [3, El E] for the 
global persistence exponent 9 g , valid for Markovian systems, from critical systems to 
phase-ordering kinetics. We have found evidence, through the breaking of the scaling 
relation ([3]) and the non-trivial form of the normalised global two-time correlator N(t, s), 
that the time-dependent global order-parameter O should in general not be described by 
a Markovian stochastic process, unless the underlying model is integrable. The example 
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of the 2D Ising model quenched to T < T c illustrates explicitly that non-markovian 
dynamics and a simple value z = 2 are independent properties of non-equilibrium critical 
systems. It is tempting to speculate that the generically non-markovian dynamics might 
be related to the fact that in phase-ordering kinetics the width of the boundaries between 
the ordered domains increases with time as w(t) ~ t K , with k — \ in the 2D Ising model 
[43] . and which might lead to some 'memory effect' in the causal interactions of the 
coarse-grained order-parameter across the interfaces. Testing this further would require 
precise simulational data for a quench right to T = (when the roughening would be 
absent) and by comparing with our results, obtained for T > 0. Further studies of the 
non-critical behaviour of the global persistence P g (t) would be very welcome. 
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